Possibility to deduce the emission time sequence 
of neutrons and protons from the neutron-proton correlation function? 
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Experimental information has been derived from the 
neutron-proton correlation function in order to deduce the 
time sequence of neutrons and protons emitted at 45° in the 
E/A = 45 MeV 58 Ni + 27 Al reaction. 
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Two-nucleon correlation functions are normally uti- 
lized to extract information on the size of the emitting 
source and on the time duration of the emission |IJ . Fur- 
thermore, when correlations of non-identical particles are 
considered, additional model-independent information on 
the emission chronology of the particles can be obtained. 

A technique to probe the emission sequence and time 
delay of ejectiles in nuclear reactions was first suggested 
for charged particle pairs Q , based on the idea that mu- 
tual Coulomb repulsion would be experienced by pairs of 
charged particles emitted with a short time delay. Com- 
parison of the velocity difference spectra with trajectory 
calculations would thus give a measure of the average par- 
ticle emission sequence ||. The technique was extended 
to any kind of interacting, non-identical particle pairs in 
the theoretical study of Ref. Q. There it was demon- 
strated that the sensitivity of the correlation function to 
the asymmetry of the distribution of the relative space- 
time coordinates of the particle emission points can be 
used to determine the differences in the mean emission 
times. This effect has been discussed for different parti- 
cle pairs, for instance pd, np irp g, K + K~ ||. The 
technique of Ref. Q has been applied in recent experi- 
mental studies where the time sequence for p and d emit- 
ted in the E/A = 50 MeV Xc + Sn reaction and for 
p, d, t, 3 He and a particles emitted in central E/A = 400 
MeV Ru + Zr collisions || has been deduced. Typical 
time delays of a few fm/c have been found in 0. 

In this paper we report on the first experimental evi- 
dence that the emission chronology of neutrons and pro- 



tons can be deduced from the neutron-proton (np) cor- 
relation function. These data come from the E/A = 45 
MeV 58 Ni + 27 Al reaction §-|l|]. After a short discussion 
about the possibility to obtain chronology information, 
clarified with the help of theoretical calculations of the 
effects, indications on the emission sequence are derived 
from the experimental ungated np correlation function as 
well as from the correlation function gated on the parallel 
velocity and on the total momentum of the nucleon pair. 

The possibility to access information from the np cor- 
relation function on the particle emission time sequence 
is connected to the different strength of the nuclear final 
state interaction experienced by different pairs. Namely, 
from a merely classical viewpoint, if there are pairs for 
which the average distance between the two particles 
(when the two particles start to interact) is smaller, these 
pairs experience a stronger interaction and exhibit an en- 
hanced correlation (or anticorrelation) strength, as com- 
pared to those pairs for which the average particle dis- 
tance is larger. Now, if there is an average time difference 
in the emission times between the two particles, there will 
also be a difference in the average distance for the two 



"classes" E„ 



> E p and E„ 



< E v 



In order to study the n—p emission chronology, one can 
calculate the np correlation function with different source 
models and parameter sets. Due to the strong final state 
interaction between n and p, the np correlation function 
shows a strong correlation for small values of the relative 
momentum q and, for some parameter sets, a weak an- 
ticorrelation for intermediate values of q. This is due to 
the fact that the correlation in the singlet channel (an- 
tiparallel spins of np) is positive while it is negative in the 
triplet channel (parallel spins of np) even if the interac- 
tion in both channels is attractive. The counterintuitive 
negative correlation in the triplet channel is due to the 
deuteron formation. Because of this process, the sample 
of np pairs is depleted and one effectively observes the 
negative correlation fll2fl . 

If a time delay in the emission of the two particles is 
introduced, the interaction is enhanced for those pairs for 
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which the average distance is smaller. This can be eas- 
ily seen if one compares the correlation function C n (q) , 
gated on pairs E n > E p , with the correlation function 
Cp(q), gated on pairs E n < E p . If the proton is emit- 
ted earlier (later) than the neutron, the ratio C n /C p will 
show a peak (dip) in the region of q where there is a 
correlation, a dip (peak) where there is an anticorrcla- 
tion, and will approach unity both for q — > (since the 
energy difference of the two emitted particles is negligi- 
ble) and q — ► oo (since modifications of the two-particle 
phase space density arising from final state interactions 
are negligible). The exact location of the peak and dip 
in the ratio depends on the source model and on the pa- 
rameter set (as does of course the correlation function 
itself). 
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FIG. 1. Calculations performed with the source model of 
Ref. [14]. Panels (a)(b): np correlation functions C(q) (solid 
dots), C n (q), C p (q) and the ratio C n /C p , from a source emit- 
ting p and n with emission time width of At = 40 fm/c 
and with the neutron emission time delayed 50 fm/c. Pan- 
els (c)(d): the same as in panels (a)(b), obtained with emis- 
sion time width of At = 80 fm/c. Panels (e)(f): the same as 
in panels (a)(b), from a "two-source model" calculation [10], 
with exponential lifetime of r p m 400±200 fm/c for protons, 
t„ ~ 600±200 fm/c for neutrons (» 30% of p and » 13% of 
n come from pre-equilibrium) . 

This effect has been calculated in Ref. [0 for differ- 
ent delays in the emission of neutrons and protons. The 
calculation for the np system assuming that protons are 
on average emitted 100 fm/c earlier than neutrons (right 
panels in Fig. 2 of Ref. [Q ) , shows C„ (q) < 1 down to q « 
15 MeV/c and a corresponding dip in the ratio C n /C p . 
The expected peak at small q, where the correlation func- 
tion shows a strong correlation, is washed out in these 
calculations, since the statistics is very low and the pa- 
rameter set generates a correlation only in the very low 



q-region (< 15 MeV/c), where the ratio C n /C p has to 
approach unity . 

To illustrate this point further, we present in Fig. 
[j] three model calculations performed with the source 
model of Ref. [H] coupled to the Koonin-Pratt formalism 
|l5| to include final state interactions. The first calcula- 
tion represents the np correlation function obtained from 
the source model of Ref. [ft4| with time width parameter 
At = 40 fm/c and with the neutron emission time delayed 
50 fm/c with respect to the proton emission time. Fig. 
|l|a shows the calculated C(q) (solid dots), C n (q) (open 
circles) and CjJq) (open squares). The ratio C n /C p is 
shown in Fig. lllb. C(q) presents a weak anticorrelation 
in the intermediate g-region between 30 and 50 MeV/c 
and a strong correlation at q < 20 MeV/c. The corre- 
lation function C n (q) gated on pairs E n > E p demon- 
strates an enhanced interaction, since the n — p average 
distance when the particles start to interact is shorter 
for these pairs. The ratio C n /C p is above unity where 
C(q) has a correlation and below unity where C(q) has 
an anticorrelation Figs. |l|c, |l|d illustrate the results of a 
similar calculation with larger width of the time distribu- 
tion At = 80 fm/c. One can see that the anticorrelation 
in C(q) is weaker and shifted to larger values of q. The 
ratio C n /C p is changed accordingly. Finally Figs. [l]e, |l|f 
present the results of the "two-source" calculation of Ref. 
JlOj l , where it is assumed that « 30% of the protons and 
« 13% of the neutrons come from pre-equilibrium emis- 
sion. The parameter set used for this calculation yields 
emission time distributions with exponential lifetime val- 
ues of t p w 400±200 fm/c for protons and r n w 600±200 
fm/c for neutrons. The ratio C n /C q remains close to 
unity for all values of q. 

The experimental np correlation function from the 
E/A = 45 MeV 58 Ni + 27 Al reaction, obtained with a 
proton detection array placed at 45° and a neutron setup 
of liquid scintillators placed at 25°, 45° and 90°, was pre- 
sented in Ref. [|l0| and is shown by the solid dots in Figs. 
§a, [2)3. Kinetic energy thresholds of 2-50 MeV for neu- 
trons and 6-50 MeV for protons are applied to the data. 
The normalization is defined in the region q = 80-120 
MeV/c. C(q) shows a clear correlation for q < 20 MeV/c 
and a very small anticorrelation for 20< q < 70 MeV/c. 

In order to study the n — p emission chronology, the 
correlation function C n (q), gated on pairs E n > E p , and 
C p (q), gated on pairs E n < E p , are presented in Fig. ^ 
together with the ratio C n /C p . E n and E p are calculated 
in the reference frame of a source moving with a velocity 
of 0.16c pi[ . A single normalization constant, calculated 
from the ungated correlation function, is utilized for both 
energy-gated correlation functions C n (q) and C p (q). 

The data show a small experimental signal. C n (q) is 
slightly enhanced with respect to C(q) for all values of 
q < 60 MeV/c (Fig. ||a). Correspondingly C p (q) is sup- 
pressed in the same g-region (Fig. 0b). The ratio C n /C p 
is above unity for q < 70 MeV/c (Fig. ||c). The signal 
has proved to be rather stable. When testing it against 
the accuracy of the energy calibrations, it turns out that 
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a statistical uncertainty of up to 20-30% does not have 
any appreciable effects on the C n /C q ratio. On the other 
hand, a systematic error of the same size might signif- 
icantly alter the final results, but we believe that the 
accuracy of our energy calibrations is within 5-10% p0[. 
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FIG. 2. Experimental ungated np correlation function 
C(q), from the E/A = 45 MeV 58 Ni + 27 Al reaction, (solid 
dots in panels (a)(b)) compared to: Panel (a), open circles: 
C n (q), constructed from pairs of type E n > E p . Panel (b), 
open squares: C p (q), constructed from pairs of type E n < E v . 
The ratio C n /C p is shown in panel (c). 

As discussed above, the enhancement of C n (q) where 
C(q) has a correlation (q < 20 MeV/c) might indicate 
that for those pairs that contribute to the region q < 20 
MeV/c protons are on average emitted earlier than neu- 
trons. The enhancement of C n (q) in the region 20 < q < 
60 MeV/c might instead indicate the opposite emission 
time sequence as C(q) has an anticorrelation there. How- 
ever, one should keep in mind that the inclusive corre- 
lation function reflects a convolution of different sources 
and origins of emission thus it is not reasonable to expect 
that the information from the inclusive data should pro- 
vide us with a well defined and unique np emission time 
sequence. 

In order to disentangle the different effects and mini- 
mize such averaging, various kinds of gates can be applied 
to the inclusive data. We shall now present correlation 
functions with additional gates and we shall see that in- 
deed, despite the poor statistics, the gated correlation 
functions yield a ratio C n /C p in closer accordance with 
the qualitative appearance of Fig. [j]. 
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FIG. 3. The same as in Fig. 2, for low-ny (< 0.14c) np 
pairs (panels (a)(b)(c)) and for high-Wii (> 0.14c) np pairs 
(panels (d)(e)(f)). 

Fig. ID illustrates how the results of Fig. || are mod- 
ified when cuts on the parallel velocity of the particles 
are applied. The effects of these cuts, that are defined 
on the basis of BUU calculations Eh and of the available 
statistics, are to enhance (suppress) the quasi-projectile 
source by selecting nucleons with large (small) values of 
the parallel velocity. Figs. [|a, ||b, ||c correspond to a low- 
i>H gate and Figs. ||d, ||e, [|f to a high-^ii gate. One can 
notice that the behavior of the correlation functions and 
of C n /Cp ratio is quite different for the two different event 
selections. "Projectile-like" events (i.e. high-Wii selection) 
exhibit an enhanced correlation function strength and an 
enhanced C n /C p ratio in correspondence to the correla- 
tion (q < 25 MeV/c). "Target-like" events (i.e. low-uii 
selection) exhibit a slight suppression in the correlation 
function strength (q < 20 MeV/c), a weak anticorrela- 
tion for q « 20-50 MeV/c and a C n /C p ratio quite close 
to unity in the correlation region. 

Gating on the momentum (or energy) of the parti- 
cle pair in the source system is particularly useful to 
investigate the time-scale of the emission and to study 
the interplay between dynamical and statistical effects in 
particle emission. The np correlation function gated on 
high-P tot pairs has already been presented in Ref. ]l0| |. 
An enhancement in the correlation function strength was 
observed indicating a smaller space-time extent of the 
pre-equilibrium source of emission p0[ . Fig. |] presents 
total-momentum-gated correlation functions. A low-P fot 
gate is applied to all results shown in Figs. ||a, [|b, |]c and 
a high-P tot gate in Figs. |d, (§3, gf. 

The effects of the gates on the correlation functions 
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and on the C n /C p ratio are similar to those observed in 
Fig. H High-total-momentum selected events exhibit an 
enhanced correlation function strength and an enhanced 
Cn/C p ratio in correspondence to the correlation (q < 25 
MeV/c). Low-total- momentum selected events exhibit 
a slight suppression in the correlation function strength 
(q < 20 MeV/c), a weak anticorrelation in the intermedi- 
ate g-region and a C n /C p ratio that does not appreciably 
deviate from unity at any values of q. 
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FIG. 4. The same as in Fig. 2, for low-P iot (< 300 MeV/c) 
np pairs (panels (a)(b)(c)) and for high- Ptot (> 300 MeV/c) 
np pairs (panels (d)(e)(f)). Ptot is calculated in the reference 
frame of a source moving with v 30urce — 0.16 c. 

The enhancement of the C n /C p ratio in the \ow-q re- 
gion, observed in Figs. ||f and Qf, could be interpreted as 
an indication that for those np pairs, protons are on the 
average emitted earlier than neutrons. 

The result of a shorter average emission time for pro- 
tons seems to corroborate our previous findings from nn 
and pp data compared to a two-component statistical 
model |l(J. However, the calculation of the C n /C p ra- 
tio with the "two-source model" parameters used in Ref. 
JTo| yields a negligible effect (shown in Fig. |l]f) as com- 
pared to the experimental data (Figs. ||f and~[|f). While 
this might indicate that the initial time distribution of 
the "two-source model" calculation is too broad and/or 
wrongly shaped, one should also keep in mind that there 
could be alternative explanations. The model of jl0| ne- 
glects any residual interaction between the emitted par- 
ticles and the reaction zone. In particular, the Coulomb 
interaction between the proton and the emitting source 
- which could increase the average distance between the 
two particles if the proton is emitted first - has not been 



taken into account pj, although, according to Ref. [[!}, 
the influence of the Coulomb field of the residual nucleus 
on the proton is expected to be minimal. Alternative ex- 
planations might be found in reaction mechanisms that 
are neglected in our calculations, such as delayed feed- 
ing from the decay of excited primary fragments. It thus 
appears as if a refined theoretical treatment would be 
welcome to confirm or not the interpretation presented 
here. 

In conclusion, the experimental results from differently 
gated correlation functions are in qualitative agreement 
with a classical argumentation and calculations using 
the Koonin-Pratt formalism JIB] , and support the in- 
terpretation that in high-parallel- velocity and high-total- 
momentum selected events, which enhance projectile-like 
and/or intermediate velocity sources, the proton is on 
average emitted earlier than the neutron. Although al- 
ternative explanations of the enhancement in the C n /C p 
ratio cannot be ruled out at this time, the experimen- 
tal evidence of such an enhancement, presented here for 
the first time, is by itself intriguing and deserves further 
experimental and theoretical investigation. 
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